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MASS TRANSFER CYCLES IN CLOSE BINARIES WITH EVOLVED 

COMPANIONS 


A. R. King, J. Frank, ^ U. Kolb ^ and H. Ritter ^ 

ABSTRACT 

We give a global analysis of mass transfer variations in low-mass X-ray binaries 
and cataclysmic variables whose evolution is driven by the nuclear expansion of the 
secondary star. We show that limit cycles caused by irradiation of the secondary by 
the accreting primary are possible in a large class of these binaries. In the high state 
the companion transfers a large fraction of its envelope mass on a thermal timescale. In 
most cases this implies super-Eddington transfer rates, and would thus probably lead to 
common-envelope evolution and the formation of an ultrashort-period binary. Observed 
systems with (sub)giant secondaries stabilize themselves against this possibility either 
by being transient, or by shielding the secondary from irradiation in some way. 


Subject headings: accretion — instabilities — stars: binaries: close — stars: cataclysmic 
variables — 


1. INTRODUCTION 

Semidetached binaries in which a compact object (white dwarf, neutron star or black hole) 
accretes material via Roche lobe overflow from a companion on or near the main sequence are of 
great interest in current astrophysics. The evolution of such systems is driven by orbital angular 
momentum losses via gravitational radiation and magnetic braking (see e.g. King 1988 for a 
review). Thus most properties of the binary, particularly the mean mass transfer rate, depend 
essentially only on the secondary mass M 2 and change on the timescale tj 10® - 10^ yr for 
angular momentum loss. In a recent paper (King, Frank, Kolb &: Ritter 1996a, henceforth Paper I) 
we discussed the conditions under which mass transfer can vary cyclically about the evolutionary 
mean in these systems by developing a general formalism allowing one to study the stability of 
mass transfer in systems driven by angular momentum losses. The existence of cycles is required to 
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account for the wide dispersion of e.g. mass transfer rates at a given orbital period. We concluded 
that the most likely cause of such cycles is weak irradiation of the companion star by the accreting 
component. The conditions for this appear to be fulfilled in a large class of cataclysmic variables 
(CVs), in which the accreting star is a white dwarf. 

Of course semidetached evolution is not restricted to systems with main-sequence companions. 
In systems with orbital periods P ^ 1 day the orbital evolution and hence the rate of mass transfer 
is either determined or strongly influenced by nuclear evolution of the companion. In this paper 
we discuss the stability of mass transfer in systems containing a giant or a subgiant companion 
and consider the possibility of irradiation-driven mass transfer cycles similar to those thought to 
exist in CVs. For this purpose we generalize the analysis of Paper I to include the effects of nuclear 
evolution on the radius variations by using a simple core-envelope model for the companion. This 
is a good representation of systems with low-mass giant secondaries, which constitute the great 
majority of long-period compact binaries. However this description does not apply to the recently- 
discovered black-hole transient GRO J1655-40, where the companion star appears to be crossing 
the Hertzsprung gap (Orosz & Bailyn, 1996). 


2. GLOBAL ANALYSIS OF MASS TRANSFER VARIATIONS 

In this section we follow closely the formalism developed in Paper I, casting the equations 
governing time-dependent mass transfer in a semidetached binary in a form suitable for global 
analysis including the effects of nuclear evolution. We restrict our analysis to the effect of variations 
in the radius i ?2 of the lobe-filling star on the transfer rate, as this is the simplest way of causing 
mass transfer variations (other ways can easily be accomodated: Paper I gives an example). Note 
that this part of the analysis is quite general, in that there is no presumption at this point that the 
radius variations result from irradiation. 

Radius changes can result from dynamical, thermal or secular processes. For example, local 
adjustments in the structure of the atmospheric layers can take place on the shortest (dynamical) 
timescale, while adjustments of the stellar radius in response to secular mass loss and nuclear 
evolution occur on the longest (secular) driving timescale tdr dehned more precisely below. We 
represent the radius variation as 

+ KttiR2,M2) + . ( 1 ) 

Here Q is the adiabatic mass-radius exponent (~ —1/3 for a fully convective star), Kth{R 2 , M 2 ) 
represents radius variations due to thermal relaxation and irradiation of the star by the primary, 
and iPnuc = (din R 2 /dt)nuc represents the secular changes due to nuclear evolution. We considered 
examples of specific forms of K{R 2 ,M 2 ) in Paper I and will discuss these again later. The change 



- 3 - 


of the mass transfer rate is given by 


M 2 


M2 , ■ ■ , ■ -R2 




( 2 ) 


where i?L is the critical Roche radius, and H is the pressure scale height in the secondary star’s 
atmosphere. The approximation given by (§) is justified since even for giant companions |i ?2 ~ 
RlI << Rl) R 2 - The response of the Roche radius to mass loss is described by 


Rl 


( 3 ) 


where Cr is a function of the mass ratio M 2 /Mi, approximately given by Cr ~ 2 M 2 /M 1 — 5/3 for 
conservative mass transfer, and J is the rate of loss of orbital angular momentum. Inserting (|l|) 
and @ into (|2|), we obtain 


M 2 




Mo 

Cr)t^ + -Rth + Rinuc 

IVI2 



( 4 ) 


For the more general case considered here the evolution of the binary is driven by the combined 
effects of nuclear evolution and angular momentum losses. We introduce the effective driving 
timescale tdr, as follows. 


1 

tdr 


1 2 

^nuc tj 


( 5 ) 


where = -Riiuc i® nuclear timescale, and tj = —J/J is the timescale for angular momentum 
losses. The nuclear timescale is a strong function of the core mass, and is ~ 10^ — 10® yr for giants 
and much longer for low mass main sequence dwarfs. The angular momentum loss time scale tj 
is typically I'N-' 10® - 10® yr for main sequence companions, but could be much longer for evolved 
companions. Thus orbital evolution is driven mainly by angular momentum losses in systems with 
main sequence companions (e.g. CVs) and nuclear processes in systems with evolved companions. 
The introduction of tdr allows us to treat both cases simultaneously. We define a dimensionless 
mass transfer rate 


-M2 

(—M2)ad 



CR)idr ) 


( 6 ) 


where (—M 2 )ad is the adiabatic mass transfer rate i.e. the steady rate implied by equation (^) with 
RTth = 0. A system undergoing stable mass transfer will typically do so at x ~ 1 (see Sect. 3). We 
also introduce the dimensionless stellar radius 





- 4 - 


r = R2/Re, (7) 

where Re = Re{M 2 , Me) is the radius of the secondary in thermal equilibrium for a given total 
mass M 2 and a core mass Me- The radii and luminosities of lower giant-branch stars are virtually 
independent of the mass of the envelope but depend strongly on the mass of the degenerate helium 
core (Refsdal &: Weigert 1970). We allow Re formally to depend on the total mass because then the 
equations obtained are identical in form to those derived in Paper I and have wider applicability. 

As we discuss in Sect. 3, the secular equilibrium radius of the companion under stable mass 
transfer differs slightly from Re- Although the equilibrium radius Re is only attained in the absence 
of mass transfer, we write formally 


Re 


C^ + K 

M2 ^ ““ ’ 


( 8 ) 


where Ce is defined by the above equation. Taking the time derivative of equation (|^, and neglecting 
secular variations of 77, we obtain 


-M2 


M 2 . M 2 X ■ d In [(Cs - CR)idr] 

(C-CR)tdr'" (C-CR)tdr ' dt 


(9) 


Finally, we introduce e = H/Re, which is typically ~ 10“^ for lower giant branch stars (see 
below). We define the dimensionless thermal relaxation function p{r,x) = ATth^dr; the dimen¬ 
sionless time r = t/(etdr)- With these assumptions and definitions (|9|) becomes 


dx 

— = rx[l + p{r, x) — xj -|- ex 
dr 


+ tdi 


d In [(Cs - CR)idr 


( 10 ) 


L(C-Cr)’^"^ dt 

Note that since the terms inside the second set of square brackets are of order unity and are 
multiplied by e, they can be safely neglected. Thus (IT) reduces to 


— = rx[l -|- p{r, x) — x] . (11) 

dr 

With the variations of the equilibrium radius given by equation (^, and using the above conventions, 
the radius equation (|I]) becomes 


dr 

dr 


= er 


p{r,x) - - 


( 12 ) 


where A = (C“CR)/(Cs“Ce)- The equations (|lli fl^ describing radius and mass transfer variations 
are identical in form to those derived in Paper I. Thus we can take over from Paper I all the results 
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for the phase plane motion of the system, the critical cnrves r = 0, x = 0, the fixed point(s) (rg, xq) 
at the intersection(s) of these critical curves, the stability anaylsis for the fixed point(s), and the 
conditions for limit cycles, with the modification that now is a more general driving time which 
inclndes both systemic angular momentum losses and nnclear evolution. In the limit tnuc —>■ oo we 
recover the case studied in Paper I. 

From equations ® 0) it is easy to see that the stationary (secular mean) mass transfer rate 
at the fixed point xq = (Cs — Cr)/ (Ce ~ Cr) depends on the properties of the companion and the 
driving rate. In particular it is independent of p{r,x) (as it must be); and in most realistic cases 
Xq ~ 1. The radins tq at the fixed point on the other hand is given implicitly in terms of p(r, x); 
again we typically have tq ~ 1 (see Paper I for farther details). 


3. BINARIES WITH SUBGIANT AND GIANT COMPANIONS 

In their stndy of the evolution of compact binaries containing a lower giant branch companion, 
Webbink, Rappaport and Savonije (1983) introduced simple analytic expressions for the luminosity 
and radins in terms of the core mass. For onr purposes it suffices to take only the first two terms 
of their approximate formulae for the radius 


and the luminosity 


ReiMc) = 12.55R0 


Le{Mc) — 33.1Lq 


M, 


0.25Mq 


M. 


5.1 


8.1 


.0.25Mo, 

in thermal equilibrinm. Using these equations we can readily estimate e for these stars 


e = ^ = 2.1 X 10"V^ 

I\p_ 


Mr. 


4.6 


( 


-1 


(13) 

(14) 

(15) 


,o.25M0y vM0y ’ 

where p is the mean molecular weight in the stellar atmosphere. From (00) the nnclear timescale 
for radial variations is tnuc = Mc/{5.1Mc), where Mr is calculated from the equilibrium luminosity 
assuming a hydrogen mass fraction of 0.7 and an energy yield of 6 x 10^® erg g~^. Thns 


tnuc = 1.0 X 10®yr 


Mr 


- 7.1 


T.25M0, 

The Kelvin-Helmholtz timescale is defined for the whole star as txH = GM 2 /{ReLe) so 

2/ Mr 


^KH = 7.6 X 10 yr 


( -^2 Y 


(16) 


(17) 


VM 0 y VO.25M0, 

The ratio of these timescales plays a major role in the discussion of thermal relaxation and stability. 
We define 


''IlUC t A t r\.^ 

Pmxc — ~ = 1.4 X 10 

IKH 


Mr 


\MqJ VO. 25 M 0 


6.1 


( 18 ) 








Since no explicit dependence of Rg on the total mass of the companion M 2 is allowed in 
equations (^, M), it appears that Ce = 0. However, both computed models of lower giant branch 
stars and analytic treatment of nearly fully convective stars following the approximations used by 
Kippenhahn and Weigert (1990) for the Hayashi line show a weak dependence on the total stellar 
mass M 2 . Taking a photospheric opacity of the form k = koP°'T^, and adopting a polytropic index 
n = 3/2 for the envelope, one can show that 

a + 3 


Ce — I 


(19) 


5.5a + 6+1.5 

With a = 1, 6 = 3, the above expression yields Ce = —0.4, while direct estimates from numerical 
stellar models indicate Ce ~ —0.2 to —0.3. We adopt the latter range in describing the reaction of 
the star to mass loss, although the relation (|T^ is adequate for computing the nuclear and thermal 
timescales. 

We now describe the effects of thermal relaxation and irradiation on the radius of the evolved 
companion. We do this here using a simple homologous model and making the same approximations 
as in King et ah 1995 and in Paper I with slight changes appropriate for giants: 

1) Homology relations can be used to describe the stellar structure, i.e. the effective temper¬ 
ature T* on the unirradiated portion of the stellar surface remains essentially constant (Hayashi 
line) and any stellar radius change changes the surface luminosity as L oa R^- In constrast L^uc 
remains unchanged, as given by equation (14). 


2) If a fraction s of the stellar surface is exposed to a uniform irradiating flux F, the luminosity 
of the star, i.e. the loss of energy per unit time from the interior, is reduced by the blocking 
luminosity 


F 

Lh = sL tanh {k — 
* 


( 20 ) 


where iT = aT^ is the unperturbed stellar flux and k a parameter which is adjusted to approximate 
the results of more realistic model calculations [e.g. those by Ritter et al. (1996a, 1996b)]. The 
above ansatz for is motivated by the facts that a) Lf^ = 0 F = 0, b) Li, ^ sL if F ^ F^,, and 
c) the transition between F = 0 and T ^ T* is smooth and monotonic. With these assumptions 
we get the thermal relaxation function including irradiation effects (cf. equation (8) of King et al. 
1995 and equation (38) of Paper I) 


pir,x) = - fp 


1 / ^ 
1 — s tanh — 
\Xc 


r — r 


( 21 ) 


Note that the quantity p = tdr/^KH introduced above differs from p = tj/^KH used in Paper I 
[see equation (^)]. The ratio / = tKH/^ce, with tee being the thermal timescale for the convective 
envelope, is a constant which depends on the internal structure (i.e. mainly on the ratio of envelope 
mass to the total mass of the secondary), and 


Xn = 


2(Cs - Cr)p ^2 Ri 


krj 


Ml Re 


( 22 ) 
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This is the dimensionless critical mass transfer rate which yields an irradiating flux F = F^,/k. 
Cycles typically occur provided that Xc is < the secular mean transfer rate xq. In (^) a is the 
orbital separation and r] a dimensionless efficiency factor relating the irradiating flux F to the mass 
transfer rate via 


rj_ GMi(-M 2 ) 
Svr Ria"^ 


(23) 


The critical mass transfer rate for CVs is Xc ~ O.lp/ki] ~ 1 — 10, where the value of k depends 
on the type of companion star assumed, while rj depends mainly on the radiation spectrum (typical 
values are 0.1 ^ p ^ 5,0.2 ^ k < 0.9, with r/ ~ a few percent). As this critical rate is of order 
the secular mean mass transfer rate xq ~ 1, irradiation can potentially drive cycles in CVs, as we 
found in paper I. For LMXBs with similar mass ratios, the much smaller neutron star or effective 
black-hole radius implies Xc ~ O.OOOlp/krj ~ 0.1 — 1. Thus we can potentially expect cycles in 
these systems also. This contrasts with the case of LMXBs with unevolved companions (cf Paper I), 
where the companion is too strongly irradiated to give cycles (critical transfer rates Xc << xq ~ 1). 
For LMXBs with evolved companions the large ratio p of driving to thermal timescales almost 
compensates the effect of the smaller primary radius to give Xc ~ 0.1 — 1. 


4. STABILITY OF MASS TRANSFER FROM AN EVOLVED COMPANION 


The simple thermal relaxation function given in the previous section allows us to discuss 
the stability of mass transfer in close binaries with evolved companions. More detailed calculations 
using the bipolytrope approximation (Ritter, Zhang &: Kolb 1996b) or integrating the full equations 
of stellar structure with appropriate changes in the boundary conditions (Hameury & Ritter 1996) 
produce results which to a first approximation can be represented by this analytic form with suitable 
A;,Xc ~ 1. In Paper I we noted that a necessary condition for limit cycles is that the fixed point 
is unstable, which in turn requires < 0 and > 1 simultaneously. These derivatives are easily 
calculable for the simple form of equation (^l|), namely 


Pr{ro,XQ) 


-fp 


1 — s tanh 




(24) 


and 

Px{ro,xo) = —rOsech^/"—. (25) 

Xc V Xc / 

For physically reasonable situations s < 0.5 and tq > 1, and thus clearly p^ < 0 as required. The 
second condition pa, > 1 for instability can be rewritten as in Paper I: 


s^sech^ 

Xc 



> 


Xq 

fpri 


(26) 


where the l.h.s. is a function which has a maximum value of ~ 0.448s at xq/xc ~ 0.78 and vanishes 
at both small and large xq/xc. As emphasised in Paper I the r.h.s. of equation (^6|) depends mainly 







on the type of companion and the driving mechanism. Clearly the large values of p for giant or 
subgiant companions makes these systems extremely vulnerable to the irradiation instability. Note 
that if s = 0 in (26), i.e. irradiation is not included, or the flux is somehow blocked (e.g. by a thick 
disk), the inequality is violated and mass transfer driven by nuclear evolution is stable as assumed 
in conventional treatments of these systems. 


Thus the simple model described above suggests that most binaries with irradiated lower giant 
branch-companions cannot transfer mass at a stable rate. More detailed models discussed later 
conhrm this result. While the linearized equations are adequate to discuss the stability of the fixed 
point, as in Paper I, we need the full non-linear evolution equations to understand the ultimate 
fate of the system. It turns out that the phase point representing the evolution of the system away 
from an unstable fixed point is trapped in its vicinity, and will settle into a limit cycle which is 
described in more detail below. 


5. PROPERTIES OF THE LIMIT CYCLE 


The simple thermal relaxation functions p{r, x) given by (^) for evolved companions and 
equation (38) of Paper I for main sequence companions yield approximate analytic expressions for 
a number of properties of the limit cycle. The two cases can be treated simultaneously by taking 


p(r, x) 


-fP 


1 — s tanh 






(27) 


with p as defined in this paper, and zz = 5 — 6orzz = —3 for main sequence and giant companions 
respectively. Some of the results given below depend on having p >> 1, which is satisfied very 
well by giants [for which p Pnuc, see (|T^)] but not always by main-sequence companions. Never¬ 
theless, since most of the analytic results described below can also be obtained for main-sequence 
companions, we shall quote them here too. 


From ([n| ) the critical curve i; = 0 is given by 1 -|- p{r, x) — x = 0. The chosen form of 
p{r,x) gives two branches of this curve with different slopes: a low branch with x « Xc and a 
high branch where x » Xc- For the low branch one can easily show using the approximation 
tanh(x/xc) ~ x/xc that 


XL 


fp(r'^ — r ^ "^) — 1 


(28) 


f par'llXc - 1 

The denominator in the above equation can be rewritten using (25) as pa, cosh^ (x/xc) — 1- The 
condition > 1 for instability ensures that this denominator is positive near xq; as xl must be 
finite the denominator must remain positive everywhere on the lower branch, implying a positive 
slope there. Since in general the radial variations are small departures from equilibrium, we can 
linearize around r = 1 to get X 2 , = 0atr = l-|- l/[/p(5 -|- i')]. In general the slope of this critical 
curve is given by 


/dx\ 


-Pr 
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where pr < 0 for physically relevant cases. If fp » 1, as in the case of evolved companions, then 
the slope of the lower branch at x = 0 is approximately (5 + v^xds. 

The upper branch is obtained by setting the tanh factor to unity and thus 

Xu ^1 + fpsd - fp{d - r~''~d . (30) 


For small departures from the equilibrium radius this reduces to x = 1 + fps — fp{5 + — 3s) (r — 1) 
showing that the upper branch has a large negative slope and is therefore stable (as s < 0.5 we 
have (5 + z/ — 3s) > 3.5 + u, which is > 0.5 even for giant companions with u = —3). At some 
intermediate value of the transfer rate xt there is a turning point at which (dr/dx)i=o = 0, or 
equivalently where px = 1- This turning point therefore satisfies the equation 


cosh^ 



sfprf 

Xc 


(31) 


where rt is the dimensionless radius at the turning point. In general rt and xt must be obtained 
by solving equation (|^ together with 1 + p{r, x) — x = 0, which can be easily done iteratively. 
Note also that as in Paper I, the axis x = 0 is formally part of the x = 0 curve, although it takes 
infinite time to reduce x to zero (see Fig. 1). From the considerations above we can see that the 
fixed point is unstable when x^ > xq > 0. 


The maximum expansion stage of the secondary is reached close to r^. We could estimate 
by asking for the radius at which xu = xu, but we prefer to use a physical argument based on the 
assumptions made in the model: the maximum expansion of the companion is reached when all the 
luminosity generated is emitted by the unilluminated side at the unperturbed effective temperature. 
This yields p{rt, xt) ~ 0 and 

r* « (1 - 1 + ^ , (32) 

v + b 


where Seff = stanh(xt/xc) ~ s. Formally if p(r,x) = 0 for non-vanishing illumination we again 
get the adiabatic mass transfer rate x = 1 and thus x = 0. This simple argument shows that main 
sequence companions expand relatively little under illumination while giants expand substantially: 
in general we have r* ~ 1 + O.lseg, 1 + 0.5sefr for these two type of companion. Thus with Sefr = 0.3 
we find from (^) maximum radii rt(ms) = 1.036, rt(g) = 1.195 for main-sequence and (sub)giant 
companions respectively {u = 5, —3). Physically this happens because upon expansion the nuclear 
luminosity of a main-sequence star decreases, whereas the luminosity generated by a giant does 
not depend on the radius. Thus assuming that the stars remain on the Hayashi line, a giant must 
expand until the unblocked area equals its original surface area, while a main sequence star finds a 
new irradiated equilibrium by expanding and simultaneously reducing its nuclear luminosity. Since 
Xf/ = 1 +p(xf/,r), the radius (^2|) {p = 0) has Xf/ = 1 with Seg = s, while the radius at which xu 
and Xu are equal differs from this only by terms ~ [f p{v + 5)]“^. We can then show from equation 
(^) that Xt is never very large, typically ~ a few times the adiabatic rate. 
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We can also make some general statements about the properties of the other critical curve r = 0. 
All physically plausible p{r, x) must be such that under no mass transfer and no illumination, the 
star remains at its equilibrium radius, implying p(l, 0) = 0. Therefore from equation (12) one sees 
that the critical curve r = 0 must go through the point r = l,x = 0. The slope of this curve at 
that point is also positive, but smaller than the slope of the a: = 0 curve at slightly larger r if a 
hxed point with xq > 0 exists. This is obvious from geometrical arguments, but can also be shown 
explicitly. 

A typical limit cycle is shown in Fig. 1, where we have stretched the r variable and exaggerated 
the separation between the two critical curves for the sake of clarity. Given that the r = 0 curve 
intersects the r-axis at r = 1, the entire cycle satisfies r > 1, i.e. the star is always somewhat 
oversized. The critical curves cross at the unstable fixed point (ro,xo). The four points labelled 
ABCD along the cycle identify the locations at which the phase point crosses critical curves. These 
naturally divide the cycle into four phases which we discuss in turn. At point A the mass transfer 
reaches its peak value Xy4(—M 2 )ad and we have maximum contact: i ?2 — Rl ~ Hln{xA/xo). We 
can estimate xa from the fact that A lies almost vertically above the point where xl = 0. From 
(|2^, 30) we have xa ~ fps. In physical units this rate is 


(-M2)r 


XA 


M2 


sM2 


iCs - CR)tdr {Cs - Cr)4 


(33) 


where tee is the timescale for thermal relaxation of the convective envelope, which can be signifi¬ 
cantly shorter than tKH- At point B the companion is very close to its maximum size, while the 
binary orbit has expanded so that i ?2 ^ Rl ~ xtRe- The degree of overfilling of the Roche lobe 
has been reduced to R 2 — Rl ^ Hln (xt/xo). The stellar expansion rate R 2 IR 2 = ^nuc 
slow compared with that of the lobe to sustain the high mass transfer rate which in turn drives 
the radius expansion. (This is inevitable since the transfer rate cannot indefinitely remain above 
the secular mean driven by nuclear expansion and angular momentum losses.) Thus the compan¬ 
ion rapidly loses contact and the mass transfer drops very sharply, while the star shrinks back 
towards its equilibrium radius. At point C minimum mass transfer is reached because the system 
is maximally detached, with R 2 — Rl ~ —{I’t — 1)/^) so that x is essentially zero in the low state. 
At this point R 2 is again very close to zero, and nothing will happen until combined effects of 
nuclear evolution and angular momentum losses bring the system close to contact again. At D the 
secondary is expanding slightly {R 2 /R 2 = Re/Re) so that r = 0. The companion radius R 2 is now 
within a few scale heights of Rl- This raises the transfer rate, making the companion expand more 
rapidly under irradiation, which in turn increases the transfer rate, so that the cycle restarts. 


We can also estimate the timescales for the different phases of the cycle and evaluate the total 
mass transferred in a cycle. We estimate the rise time by assuming that the mass transfer initially 
increases because the thermal imbalance due to irradiation forces the secondary into deeper contact. 
The characteristic rate of radial expansion when a fraction s is fully blocked is dlni22/dt ~ s/tce- 
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Therefore the rise time is approximately the time required to expand by In [xa/xq) scale heights 

\n{xA/xo)H tee ln(xA/xo) 

toA ~ - or tda ~- , (34) 

R 2 S XA 

where we have used our estimate of xa = sfp to obtain the final expression. During the time Iab 
spent on the high branch, the rate of expansion of the secondary decreases monotonically from the 
maximum rate estimated above, which yields the peak mass transfer rate xa-, until dlnii 2 /dt = 
1/tnuc- At that point (B) we have r = 0 and the expansion rate falls below the driving rate so that 
the system detaches, going rapidly into a low state. We adopt the following ansatz for the radius 
as a function of time 

R 2 (t) = RA + iRB-RA){l- ) . (35) 

Equating dlnii2(0)/dt = s/tce and dlnR 2 (tAB)/(R = 1/tnuc, we obtain both the characteristic 
radial expansion timescale T+ and the time tAB spent on the high branch 


in - 1)^ 


n -1 


in - 1) 


^ce In ixA/n) 


u + b 


In (xA/n) 


n -1 


\n{xA/rt) , 


where we have taken Ra ~ Re and Rb ^ XtRe- As soon as the system detaches, irradiation ceases, 
and the companion contracts rapidly while the orbit and Roche lobe remain at the size attained at 
the end of the high state. The contraction is even more rapid than the expansion because the star 
is more luminous. One can show that the characteristic thermal contraction time scale is 


T_ ^ (1 - 


leading to 


rt-l tee 

rf Seflf 




in - i)t 


ce D '®eff^dr 


^eff ^ce ^ce 

Clearly this effect is more pronounced in giant companions than in main-sequence secondaries, 
for reasons that have already been mentioned. In a few times r_ (i.e. typically a time ^ 0.5tce) 
the mass transfer formally reaches a minimum value xc ~ 0 and a long detached (low) state now 
follows while the driving tries to bring the system back to contact. The time spent in the low state 
is thus dominated by the time ten 


tcD ~ in - l)tdr ~ -TT^dr 

U + O 


(n -1) 


The total mass transferred during a cycle can now be estimated as AM 2 
which yields the simple - with hindsight perhaps obvious - result 


i — M2)adirt — l)idr — - ■^M 2 


(Cs — Cfi)(^ + 5) 


XAi-M2)s.dT+, 


( 41 ) 
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Hence the total cycle time tcyci — AM 2 /(—M 2 )o is 

rt — 1 (r^ — 1) . . 

tcycl ^ - tdr or Tcycl ~ - • (42) 

Xo exo 

Thus irradiated main-sequence stars (i^ = 5 — 6) transfer at most a few percent of their mass per 
cycle. In contrast giants {v = —3) transfer a significant fraction ~ s of their total mass M 2 in the 
high state, which may amount to most of the convective envelope. Since tee is also much shorter 
for giants we expect that if the irradiation instability is allowed to grow unchecked in such systems 
it will produce accretion rates which are super-Eddington in LMXBs. The resulting common- 
envelope evolution would probably lead to the formation of an ultrashort-period binary. Clearly 
the irradiation instability must be quenched in observed LMXBs with evolved companions. In 
Section 7 we discuss ways in which this can happen. For reference we collect together the analytic 
expressions for the various properties of the limit cycle in Table 1. 


6. NUMERICAL RESULTS 


We have integrated the evolution equations do^iii) from arbitrary initial states for a variety of 
parameter choices. In this section we present a few examples and compare them with the analytic 
estimates given in the previous section. The values quoted in the text are those obtained by 
numerical integration, with the corresponding analytic estimate in brackets. The analytic estimates 
for peak rates, rise times, and maximum expansion radius are in all cases calculated in very good 
agreement with numerical results. In general we integrated the equations for several cycles (3-6) 
and noted that the first outburst is usually somewhat untypical (slightly higher for the same radii 
and a higher peak rate) because of initial conditions, whereas later outbursts are virtually identical. 
We therefore quote values taken from later outbursts. 


Fig. 2 shows a case with core mass Me = O.ISMq, the lowest reasonable value for which 
the approximations used in Section 3 are still valid. We have also taken s = 0.5, which is perhaps 
unrealistically large unless there is significant scattering from a disc corona or similar structure. The 
luminosity and radius of the companion are at the lower end of the subgiant range and therefore the 
instability is weakest. Nevertheless a very large amplitude outburst results, with a peak transfer rate 
XA = 152.8 (154.9) which is attained rapidly, in t^a = 0.032 (0.033). The companion continues 
to expand in the high state until the turning point is reached at rt = 1.41421 (1.41421). The 
duration of the high state can be read off the graph directly and is tab = 8.1 (12.5). While the 
behaviour of r as shown in Fig. 2 appears to obey equation (^), closer examination shows that 
the expansion phase is actually faster than exponential while the contraction phase is slower than 
the assumed exponential. Nevertheless we have estimated the characteristic radial expansion and 
contraction timescales graphically obtaining r_|_ = 2.4 (2.67) and r_ = 1.0 (0.95). Despite the fact 
that equation (^) holds only approximately, the expressions based on it are better than order of 
magnitude estimates. For the case shown, the total duration of the cycle is 415.9 and the duration 
of the low state is tcd ~ 407.8 (414.2). 
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Figure 3 shows another case with a higher core mass = O.25M0 and s = 0.3, appropriate 
for illumination by a point source. The instability in this case is more violent, with shorter rise 
times To A = 0.0023 (0.0024) and higher peak rates xa = 3326 (3609), but a smaller turning radius 
rt = 1.19529 (1.19523). The total duration tab = 0.45 (0.43) of the high state is relatively short. 
The characteristic e-folding times for radial expansion r+ = 0.051 (0.054) and contraction r_ = 0.04 
(0.031) are also relatively shorter because of the higher luminosity. The analytic estimates for the 
properties of the radial variations during the limit cycle are even more accurate than in the case in 
Fig.2, because the neglected terms ~ ifps)~^ are still smaller here. The total duration of the cycle 
is 188.6 and the low state lasts for tcd ~ 188.1 (195.2). 


7. THE IRRADIATION INSTABILITY IN LOW MASS BINARIES 


We can now discuss the application of the theory developed here to various types of close 
binary encountered in nature. The only restriction is that the companions must have a significant 
convective envelope and thus a low mass M 2 ^ 1.5Mq. The instability criterion implying cycles 
given in (|^) can be rewritten as 


fps > —E — cosh^ 
ro xo 


(43) 


Here the factor xq/tq ~ 1 does not vary much, whereas fps is mainly sensitive to the type of 
companion and the mechanism driving the binary evolution, while ^ = x^jxc depends on both 
the accretor and donor type. In Fig. 4 we plot x^jxc along the abscissa (the compact object 
axis) and fps along the ordinate (the companion axis). Changing the type of companion causes 
displacements along both axes, whereas changing the primary causes only horizontal displacements. 
The stability/instability boundary plotted is simply fps = .^“^cosh^^. The locations of various 
possible evolutionary sequences for CVs and LMXBs are also shown on Fig. 4. Changing s by 
screening or scattering causes purely vertical displacements since | is not affected: clearly for any 
binary there is a value of s below which mass transfer is stable. Within the limitations of the 
simple theory developed here (i.e. the form chosen for p{r,x)) the results displayed on Fig.4 show 
the following: 


1. CVs above the period gap and above a certain period (or companion mass) can be unstable, 
depending on the value of p (see Paper I). CVs below the gap are stable. 


2. Long period CVs with companions having small core masses are stable, whereas companions 
with larger core masses likely to have even longer orbital periods are unstable. A detailed analysis 
shows that GK Per and V1017 Sgr are unstable if r/ > 0.08 and 0.04 respectively. 


3. Short-period LMXBs with main sequence or partially evolved companions are stable because 
they are so strongly irradiated that they have reached saturation {xq >> Xc)- Thus variations in x 
do not cause radius variations, eliminating the feedback necessary for instability. 
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4. Long period LMXBs with (sub)giant companions are unstable: the larger the core mass and the 
smaller the total companion mass, the more violent the instability. 


The irradiation instability may well cause the formation of ultrashort-period systems (e.g. 
the AM CVn’s) from long-period CVs and LMXBs. However, since the rise to the high state is 
so rapid it is extremely improbable that we currently observe any long-period LMXB undergoing 
cycles. We must therefore consider ways of quenching the instability in these systems. The most 
obvious possibility is screening of the companion from the irradiation, which is the basic cause 
of the instability. Screening may well result from the extensive disc coronae inferred in LMXBs 
(e.g. White &: Mason, 1985). Formally we can consider this possibility by decreasing s in our 
simulations. As can be seen from equation (^) the effect of a small s is to lower the high branch 
and thus reduce the amplitude of the cycle in x. Since the lower branch remains unaffected, the 
turning radius and hence the radial amplitude will also be reduced. This suggests that screening 
could reduce the amplitude of the cycle and increase the frequency of outbursts, because it will 
take less time to drive the system into contact again once it has detached following an outburst. 
However this argument applies only to a fixed s. If s(a:) is itself a function of the instantaneous mass 
transfer rate, as a result of a varying geometrical thickness of the accretion disk or varying optical 
depth through the (out)flow, a more gradual transition occurs. For example, if screening results in a 
rapid reduction of s beyond some critical Xgcr, the amplitude of the cycle in x is reduced, saturating 
somewhere close to Xscr, whereas the radial amplitude remains large; thus the outburst does not 
recur any sooner than in the unscreened case. We have performed some numerical experiments to 
simulate these effects and verify the above statements. These simulations and the discussion above 
imply that if Xgcr > xt, the radial amplitude of the cycle and recurrence time are not significantly 
different from the unscreened case. However, if Xscr < xt, then both the radial amplitude and the 
duration of the low state are reduced. Also the duration of the high state and the amount of mass 
transferred per cycle are decreased. Clearly if Xgcr is reduced to ~ xq we will find that the cycles 
disappear altogether. Thus efficient screening can either eliminate the cycles entirely or render 
them relatively harmless as far as the binary evolution is concerned. 


While screening must play a role in stabilizing some systems against the irradiation instability, 
a second way of quenching the instability appears to be more common. This mechanism uses the 
fact that the companion swells only when its own intrinsic luminosity is blocked by the irradiating 
flux, and that the blocking effect saturates once the latter is comparable with the intrinsic flux. 
A variable accretion rate, as seen in e.g. soft X-ray transients, severely reduces the efficiency of 
irradiation in expanding the companion: the very high irradiating flux during outbursts has no 
more effect than a much weaker value, while the star can cool between outbursts. We might thus 
expect a highly modulated accretion rate with a duty cycle d << 1 to mimic irradiation by a 
steady accretion rate a factor d smaller. This expectation is largely fulfilled, as the following simple 
calculation shows. 


We assume the dimensionless accretion rate to vary periodically in time over t^ec as 


Xacc — Xfi 


0 < t <th 


= Xl, 


ih ^ t — ^rec 


( 44 ) 
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Mass conservation requires that the dimensionless transfer rate obeys 


X = dxh + {I - d)xi = Xh 


d + 


l-d- 

^4“ 


(45) 


where d = th/Uec is the duty cycle and A = Xh/xi the amplitude of the accretion rate variation. 
The reaction of the companion to this intermittent irradiation is governed by the thermal relaxation 
function (p7|), which now becomes 


and 


P = - fp 


P = - fp 


1 — s tanh 


1 — s tanh 


(■ 


Xc d{A — 1) + ly 
c 1 


^3 _ (i^+2) 


0 <t <th 


^3 _ ^-{u+2) 


\Xc d{A — 1) + ly 

Now assuming that the cycle time tree is << tee, we can dehne a mean value 


th<t<tr 


(46) 

(47) 


< p >= 


1 



r)dt . 


(48) 


Performing the integration, we can compute the derivative < p >x which governs the stability of 
the fixed point (xo,ro). If Xh >> xi, i.e. A —>■ oo, as is characteristic of soft X~ray transients and 
dwarf novae, we have 

<p{xo,ro) >x= •^^^^° sech^ 

Xc \dxcj 

and the criterion for instability < p(xo,ro >x> 1 can be written as 


fps > 


Xo 




(50) 


This criterion is very similar to the earlier one (^) assuming steady accretion, to which it of 
course reduces as d —> 1. For xo/dxc ^ 1 the two criteria are identical. Thus on Fig. 4 the 
stability/instability boundary for d < 1 is given simply by sliding the d = 1 curve parallel to itself 
along the asymptote for small xojxc, by a displacement of — logd in xq/xc- We may thus draw a 
further conclusion from Figure 4: 


5. Typical soft X-ray transient duty cycles d < 10 ^ are probably enough to stabilize most LMXBs 
with evolved companions against the irradiation instability. 


By contrast, extremely short duty cycles d ^ 10“^ would be required for dwarf nova outbursts 
to stabilize CVs with evolved secondaries. We note that most LMXBs with periods ^ 1 d are 
transient (King, Kolb Burderi, 1996), and indeed this holds for a very large fraction of long- 
period systems also (King, Kolb & Burderi, 1996; King et ah, 1996b). We shall consider the 
application of the stability criteria to individual systems in a future paper. 
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8. CONCLUSIONS 

We have shown that irradiation of an evolved low-mass companion in LMXBs and CVs can 
drive mass transfer cycles. These cycles do not need the intervention of any further effect, unlike 
the case of main-sequence companions (Paper I), where a modest increase in the driving angular 
momentum loss rate is required in the high state if cycles are to occur in many cases. The cycles 
with evolved companions are also considerably more violent than in the main-sequence case. This 
is a direct consequence of two facts. First, the nuclear luminosity of an evolved star is insensitive 
to the stellar radius, so that blocking of the intrinsic stellar flux by irradiation requires the star to 
expand so as to maintain the same unblocked area. This leads to much larger expansions than in 
the main-sequence case, where the nuclear luminosity drops sharply as the star expands. Second, 
the ratio tdr/Ue = fp of driving to thermal timescales is much larger in evolved stars than on the 
main sequence, making the expansion very rapid. In the high state an evolved star would lose a 
significant fraction of its total envelope mass on a thermal timescale. 

In a CV with a low core mass, the implied accretion rate would probably turn the system into a 
supersoft X-ray source. In CVs with higher core masses, the nuclear burning causes the white dwarf 
to develop an extensive envelope, while in long-period LMXBs, the high state accretion rates greatly 
exceed the Eddington limit. If the instability is not quenched all these systems would undergo a 
common-envelope phase. They may merge entirely, or reappear as ultrashort-period systems like 
AM CVn (for white dwarf primaries), or helium-star LMXBs, or detached systems with low-mass 
white dwarf companions. However, there are at least two ways in which the instability can be 
quenched and the systems (particularly LMXBs) stabilized; shielding of the companion by e.g. 
an extensive accretion disc corona, and intermittent accretion. The typical duty cycles d ^ 10“^ 
observed in soft X-ray transients are short enough to stabilize them, while observed dwarf nova 
duty cycles are unable to stabilize CVs with evolved companions. 

Although both means of stabilization seem to occur in nature, it is clear that the irradiation 
instability is so violent that it must play a major role in any discussion of the evolution of CVs and 
LMXBs with evolved companions; the systems must somehow stave it off or evolve catastrophically. 
We shall consider some of the observational consequences in a future paper. 
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Table 1; Characteristic Properties of the Limit Cycle 
Parameters 


tdr = (1/tnuc + 
tee = 

Cs = {dlnR2/dlnM2)s 
Ce = (91ni?2/91nM2)e 
(r = {d In Ri^/din M 2 ) 
e = H/R2,e 

{ 5 — 6 MS donors 
—3 giant donors 
Seff = s tanh(xt/xo) ~ s ~ 0.3 


driving timescale 

thermal timescale of the convective envelope 
adiabatic mass radius exponent 
thermal equilibrium mass radius exponent 
Roche lobe mass radius exponent 

photospheric scale height in units of equlibrium radius 
(— 1 / — 3) is the radius exponent of the nuclear luminosity 

effective blocked surface fraction 


Characteristic Properties in the limit Seff^dr/^ce 1 
in physical units in dimensionless form 


(-M2)o = 
{-M2)a ~ 


M 2 1 


Ce Cr ^dr 
M 2 Seff 


Cs Cr tee 
R2,t = R2,B ~ (1 - 

AM 2 W 7 ^2 

Cs Cr 

, 1 n ^eff^dr 

^AB ~ ^ce ■ 


^eff 


ntc 


, 1 ice n 

tBC - — 3 -In —:- 

Yf 5eff tcQ 


tcD ~ (n - 1) tdr 

^tce 1 'Seff^dr 


toA 


In ■ 




Tf — I 

tcycl ~ ^dr 

Xo 


Xq = 


XA 


Cs — Cr 
C e - Cr 

Seff^dr 

tee 


Seff) 

AM 2 ^ Tt- I 
M 2 Cs - Cr 

n - 1 , XA 


tab 


In — 

e XA n 


rt-1 

tbc ^ -3- 

r^exA 

n-I 


InYxA 


TCD 

tda 

Rrycl 


1 , XA 
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Fig. 1.— Phase plane for the ODE system (H, ^). The two critical curves where x = 0 (dot-dot-dashed) 
and f = 0 (dashed) are shown for a typical p{r,x). These curves intersect at the fixed point (ro,xo) and 
divide the (r, x) plane into four regions in which the motion of the system point is indicated by the arrows. 
The limit cycle intersects the critical curves at the points ABCD giving the four phases of the cycle discussed 
in the text: a high state AB during which the companion expands, a moderately rapid contraction phase 
BC, a long low state CD, and an extremely fast rise DA to peak mass transfer. 
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Fig. 2.— The dimensionless mass transfer x and radius r evolution during a typical outburst for a subgiant 
companion with core mass Me = O.ISM© and a large illuminated fraction s = 0.5. The calculations assumed 
a fixed tj = 10^° yr, a neutron star primary with Mi = 1.4M0 and i?i = 10® cm, M 2 = O.SMq, e = 0.001, 
kr] = 0.01, C = -1/3, Ce = -0.2, and Cr = 2 M 2 /M 1 - 5/3. 
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Fig. 3.— The dimensionless mass transfer x and radius r evolution during a typical outburst for a giant 
companion with core mass Me = O.25M0 and an illuminated fraction s = 0.3, consistent with a point source. 
The calculations assumed a fixed tj = 10^° yr, a neutron star primary with Mi = l.dM© and i?i = 10® cm, 
Ma = O.5M0, e = 0.001, kr] = 0.01, ^ = -1/3, Ce = -0.2, and Cr = 2 M 2 /M 1 - 5/3. 
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Fig. 4.— The approximate stability boundary given by equation (M) in the (xq/xc, /ps)-plane, for various 
values of the accretion duty cycle d. Representative evolutionary sequences for different types of low-mass 
binaries are also indicated. These include black holes, neutron stars and white dwarfs with main-sequence 
and low-mass giant companions (denoted BLMXB, NLMXB, CV respectively). We took black-hole, neutron 
star and white dwarf masses Mi = IOMq, 1.4M0 and I.OMq, an irradiated fractional area s = 0.3, krj = 0.1 
in short-period systems, and krj = 0.01 in long-period systems. The latter were evolved from an initial core 
mass Mci = O.ISMq and total secondary mass M 2 i = O.SMq. We see that steady LMXBs with evolved 
companions are prone to the irradiation instability, while CVs with giant companions are unstable if their 
core masses are sufficiently large (see text). Typical soft X-ray transient duty cycles d < 10“^ quench the 
irradiation instability in LMXBs, while typical dwarf nova duty cycles cannot do this for CVs. We stress 
that while the sequences shown here are representative, the stability or otherwise of any given individual 





